Infrared optical absorption spectra of CuO single crystals: Fermion-spinon band and 
dimensional crossover of the antiferromagnetic order 
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We have obtained mid-infrared optical absorption spectra of the 5 = 1/2 quasi one-dimensional 
CuO using polarized transmission measurement and interpreted the spectra in terms of phonon 
assisted magnetic excitations. When the electric field is parallel to the main antiferromagnetic 
direction a A shaped peak is observed with the maximum at u = 0.23 eV which is attributed 
to spinons along Cu-O chains. At low temperatures in the antiferromagnetic phase another peak 
appears at uj — 0.16 eV which is attributed to two-magnon absorption but the spinon peak remains. 
This behavior is interpreted as due to a dimensional crossover where the low temperature three- 
dimensional magnetic phase keeps short range characteristics of a one-dimensional magnet. 
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Properties of one dimensional systems are radically dif- 
ferent from those of higher dimension. For Heisenberg 
one-dimensional (ID) antiferromagnets (AF) the ground 
state is a spin liquid with fermionic excitations called 
spinons rather than the usual, propagating spin waves, 
of ordered higher dimensional systems. A fundamental 
question in condensed matter physics is what of these 
anomalous properties and to which extent can survive in 
anisotropic but higher dimensional systems. For a quasi- 
1D system formed by chains bridged by small perpen- 
dicular interactions the ground state is expected to be 
ordered, however a spin liquid may still be a good approx- 
imation for the ground state. Furthermore 2D quantum 
antiferromagnets, although ordered at T = 0, are be- 
lieved to be close to a spin liquid ground state. This 
has been widely studied in connection with high tem- 
perature superconductors as it has been proposed that 
superconductivity emerges from a spin liquid 1 -. Even in 
an ordered state, intermediate to short range spin liquid 
characteristics may reflect on the response at finite fre- 
quencies, as proposed by Anderson and collaborators to 
explain puzzling side bands in optical spectra of 2D AF 
cuprates 2 -. 

In this communication we study the mid-infrared spec- 
trum of CuO (tenorite) due to magnetic excitations 
infrared (IR) active by the Lorenzana and Sawatzky 
(LS) phonon-assisted mechanism 3 ' 4 . CuO is a quasi- 
1D spin-1/2 compound (Cu 2+ ) which becomes a 3D 
antiferromagnets at low temperatures T < Tjvi = 
210K. It has an intermediate phase at T/vi < T < Tjv2 = 
230K which is a spiral and has attracted recent atten- 
tion for being a high-T c multiferroio^. We show that the 
mid-IR spectrum of CuO at T > Tjv2 can be interpreted 
as due to spinon excitations as found for S^CuOs, the 



well known ID S=l/2 AF material^. In the ordered 
phase (T < T/V2) a feature due to the propagation of 
two spin waves appears similar to the bimagnon peak of 
2D cuprates 3,4 . The high energy side bands however re- 
main strong and evolve continuously from the fermionic 
spectrum in strong analogy with the scenario of Ref. 2 in 
which fermionic excitations persist in an ordered system. 

The tenorite CuO single crystal with the monoclinic 
lattice structure was grown using the floating zone 
method^. At T < T Nl the Cu 2+ spins (S=l/2) are paral- 
lel to the [010] direction (=6-axis) pointing normal to the 
ac-plane (see inset of Fig. Q] (b)). Within the plane they 
show antiferromagnetic alignment along the [101] direc- 
tion and ferromagnetic order along the [101] direction. 
We label the two magnetic directions as A=[10l] and 
F=[101] respectively. To probe IR response of the three 
principle axes A,&, and F, we have prepared three samples 
with (A,F),(A,&) and (6,F) faces respectively with typical 
area of 5mm x 5mm. For IR transmission measurement, 
they were mechanically polished to ~50 /Ltm thickness. 
Incident IR radiation was linearly polarized such that 
the electric field E was parallel to the magnetic axes. 
For example in the (A,F) sample the transmission T(w) 
was measured with the IR electric field E and magnetic 
field B set as (E//A, B//F) and (E//F, B//A). With 
the three samples we collect the full set of data, (E, B)// 
(A,F), (A,6),(F,6),(F,A),(6,F) and (6,A). From them we 
find that T(uj) is determined by E but not dependent on 
B. For example, T(ui)'s for (E, B)// (A,F) and (A,6) 
were identical. It holds when E is polarized along C and 
b axes as well. In the monoclinic CuO, A and F are not 
orthogonal but form an angle of 99.54°. Therefore it is 
not possible to get a "pure" F or A response but there is 
"leakage" in T(cu) although the effect is small. The opti- 
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FIG. 1: (Color online) (a) Mid-infrared conductivity of 
CuO crystal when IR electric field is polarized along the 
A-axis (E//A), F-axis (E//F) and &-axis (E//6) respec- 
tively, (b) Temperature dependent optical conductivity <xi(w) 
for E//A mode. The inset shows the spin ordering at T< Tjvj 
in the ac-plane. Here A and F represent the antiferromagnetic 
(A) and ferromagnetic (F) axes, respectively. 



cal conductivity ai(u>) for the three modes E//A, E/ /F, 
and E//6 were extracted from T(ui) using the standard 
Fresnel formula analysis. 

Fig. Ufa) shows the optical conductivity of CuO taken 
at T = 80K when E is parallel to A, F, and b axes. In 
E//A mode, the mid-IR absorption consists of a broad 
A shaped peak labeled II with a much narrow peak su- 
perposed on the low energy part labeled I. These Peak I 
and II are absent in the ferromagnetic axis E//F. This 
is consistent with an electric dipole operator— containing 
a term Sj.Sj which excites AF aligned spins in sites 
In the E//b polarization we find an absorption bands 
with an intermediate strength. The rapid rise of o\ (cj) 
at w < 140 meV is due to the optical phonons. The 
infrared phonons of CuO crystal were studied by Kuz- 
menko et al. from reflectivity measurement^ according 
to which the highest energy phonon peak lies at 65 meV 
for E//A. We have checked that the tail of such phonon 
overlaps the rise of <7i(w) below 140 meV. The strength 
of the phonon at its center frequency is about 800(0- 
cm) -1 , much greater than that of the mid-IR peak. It 
shows that although both features are E-active, the mid- 
IR absorption has a different activation origin from the 
phonon. 

According to the LS mechanism magnetic excitations 
which are not directly dipole active, become visible in the 
IR spectra thanks to the assistance of an optical phonon. 
At lowest order two magnons (or two spinons in ID) plus 
a phonon are created by the absorbed photon. Generally 
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FIG. 2: (Color online) The left panel shows the spinon con- 
tinuum (shaded region) in the u-q plane. The right panel 
shows the line shape that results after momentum integration 
(neglecting the phonon). 



the spectrum is similar to the magnetic Raman spectrum 
with the difference that because the phonon carries mo- 
mentum, an extra integration over the total momentum 
of the magnon (or spinon) pair, q is required, instead of 
the pair having zero total momentum as in Raman. Thus 
the spectrum may have sharp structures due to Van Hove 
singularities which provide precious information on the 
excitations in the system. Another important difference 
with Raman is that the spectrum at T — is shifted 
by the phonon frequency thus appearing as a phonon 
side band. The relevant phonon is usually a stretching 
mode phonon whose frequency can be estimated from 
the highest optical phonon. In the following we take 
ujph = 65 meV 10 for the E//A mode which we analyze 
in detail. The AF chains can also contribute in the b 
direction involving a b polarized buckling phonon which 
can add to interchain magnetic contributions involving 
other phonons. Thus the E//6 line shape has a mixed 
intrachain -interchain character. This more complicated 
case probably relevant for the multiferroic effect will be 
analyzed elsewhere. 

We note that peak I, in E//A mode, is similar to the 
two-magnon peak observed in 2D cuprates at ~ OAeV 
which is dominated by bimagnon o 3 ' 4 ' 11 ^ 2 while peak II 
resembles the A shaped peak observed in ID cuprates^ 
and well understood in terms of spinon absorption 9 . In 
a purely ID AF Heisenberg system with exchange inter- 
action J, the dipole operator connects the ground state 
with the two spinon continuum which is limited by the 
des Cloizeaux-Pearson dispersion relation^ 3 - (the shaded 
region in the left panel of Fig. [2|) . At T = the spectrum 
of spinon excitations can be computed virtually exactly 9 
while at finite temperatures numerical and analytical re- 
sults are available^. The lower edge of the continuum 
has a saddle at momentum q = n/2 (setting the lat- 
tice spacing a = 1) with energy uj(q = n/2) = TtJ/2. 
When integrated over momentum at T — a Van Hove 
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FIG. 3: (Color online) Full lines: Temperature dependent 
changes of experimental (full curves) and theoretical (dashed 
lines) ai(u>) of the E//A polarization for (from top to bot- 
tom) T = 80K, 150K, 200K, 220K, 260K and 295K. Curves 
where shifted by O.lf^cm -1 for clarity. A common constant 
background was subtracted to the experimental data. For the 
theory we used an effective J varying linearly with tempera- 
ture with J(295K) = 94 meV and J(80K) = 100 meV and 
with similar small adjustments to the intensity scale. 

singularity appears at u> = Lu p h + 7rJ/2 corresponding 
to the maximum of the absorption^ (c.f. Fig. [21 right 
panel). Increasing the temperature the Van Hove sin- 
gularity broadens^ but the maximum remains close by. 

Fig. [3] compares the experimental data with the two 
spinon theoretical ansatz of Ref,-^ for a purely ID sys- 
tem at different temperatures. Surprisingly despite the 
substantial 3D character of the magnetism in CuO, the 
main feature is in good agreement with the theory with 
a weakly varying temperature dependent effective J. At 
high temperatures the main difference is the high energy 
tail which may be attributed to higher multispinon or 
multiphonon/multispinon processes. Thus above the or- 
dering temperatures the line shape is well explained by 
the fermionic spinon excitations of a pure ID system. In- 
deed at high enough temperature spins on different chains 
are uncorrelated and the small interchain interactions av- 
erage to zero. As the temperature is lowered the effect of 
the interchain interactions becomes to be relevant. One 
effect is to renormalize upwards the effective J along the 
chains due to indirect exchange processes. This explains 
the weak temperature dependence of J needed to fit the 
experiment. Another effect will be to drive 3D magnetic 
order as discussed below. 

In order to understand the origin of peak I we plot in 
Fig. Q] the change of optical conductivity with respect to 
the room-T data. Note that peak I and peak II show 
different T-dependences; while peak II evolves gradually 
in this temperature range, consistently with the fact that 
KbT « J, peak I develops below Tjvi trough a reor- 
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FIG. 4: (Color online) Optical conductivity difference 
Acri (T) = (7i (T) - (Ti (295 K) . 

ganization of the low energy part of the spinon contin- 
uum. This can be seen more clearly in Fig.[5ja) where we 
plot the change in oscillator strength of the two peaks as 
AS(T) — J A<Ji(u>)dui where the integral is taken for the 
relevant ranges 153 meV ~ 186 meV for I and 186 meV 
~ 260 meV for II, respectively. Peak I grows reorganiz- 
ing the spinon spectral weight which is transferred from 
low energies to high energies. In addition the whole mag- 
netic line shape gains weight as the temperature is low- 
ered. This is due to transfer of spectral weight from the 
phonon region to the magnetic region as the AF ordering 
builds up and the excitation of antialigned spins become 
more likely. Notice that these two processes keep the 
total sum rule practically constant [c.f. Fig. [5{c)]. On 
passing Tjvi we notice that an anomaly appears on the 
temperature dependence of peak II when the spiral de- 
velops at Tn2, which may be related to the mutiferroic 
effect. 

The temperature dependence of peak I shows that it re- 
quires long-range order where spin waves became well de- 
fined propagating excitations. The dipole operator pro- 
duces two spin waves in nearest neighbor sites which in- 
teract strongly^. A rough estimate of the position of 
the peak can be obtain from the energetic cost to flip 
two neighboring spin, along the AF chain, in the Ising 
limit. Taking into account only J along the chains one 
expects the peak to appear at u p h + J = 0.159 eV in 
good agrement with the observation. Clearly the smaller 
interactions which drive the 3D AF order can not be ne- 
glected and at first sight they would spoil the agreement. 
Their effect would be to increase (decrease) the estimate 
providing they are non-frustrating (frustrating) with the 
non-frustrating effect expecting to dominate. To some 
extent this is already taken into account by the renor- 
malization of J however some hardening is still expected. 
On the other hand more realistic interacting spin wave 
theory (ISWT) computations give generally a smaller re- 
sult than the Ising estimate so a cancellation of errors 
may explain the surprisingly good agreement of the sim- 
ple estimate uj p h + J- For example in 2D one obtains 
AEjswt = 2.7J while AE Ismg = 3J. A precise ISWT 
computation in the present case is not possible because 



4 




3 232 
30 

% 20 
T 10 

a „ 

T 

o 

-10 

<3 -20 
-30 





■ aS(coxd c ) 


(c): 

-\ I 1 


-H 1 1 1 I— 

- aS(co<co c ) 


a) c = 153meV - 



100 



150 200 
T(K) 



250 300 



FIG. 5: (Color online) (a) Optical spectral weight change 
of the peak I and peak II estimated from AS(T) — 
j Aai(T, uS)dui. (b) Temperature dependence of the max- 
imum of peak II. (c) AS(T) for the low frequency region 
80meV < uj < 153meV and for the high frequency mid-IR 
peaks at 153meV < uj < 375meV. 



the structure is much more complicated and the inter- 
chain interactions are not precisely known. 

In the ordered temperature range Acti(w) exhibits two 
weak sideband peaks at higher frequencies as indicated 
by the arrows. Their peak energies, 0.276eV and 0.334eV 
are close to ui p h + 2 J = 0.216eV and w p j, + 3 J = 0.316eV 
consistent with multi-spin-wave excitations as observed 
in 2D cuprates. Also the position of peak II has an 
anomaly at the ordering temperature and hardens below 
Tjvi. All this shows that the reorganization of the spinon 
spectral weight due to the emergence of the spin waves 
occurs also at high energies. Spinons and spin- waves mix 
but conserve to a surprising extent their own character. 

The value found for J ~ 95 — 100 meV is larger 
than a previous "rough" neutron scattering estimate 5 
(67±20 meV), a susceptibility measurement^ (73± 
3 meV) and self-interaction-free LDA (52 meV)^ but 



smaller than a recent LDA study (126 meV or 118 meV 
depending on detail of the method)^!. Notice also that 
our measurements weight more heavily the high energy 
part of the spectrum which may be affected differently by 
the smaller interactions than the previous experiments. 

To conclude we have shown that CuO can be described 
at high temperatures as a ID quantum antiferromagnet 
with fermionic excitations. As the temperature is low- 
ered a dimensional crossover occurs and the system be- 
comes ordered with well defined propagating spin-waves 
but with a strong remnant of the spinon band. A sim- 
ilar phenomena has been seen in KCUF3 using neutron 
scattering where the spin waves appeared as a very subtle 
reorganization of the low energy spinon spectral weight 18 . 
Here the effect is much more dramatic challenging theo- 
ries of the crossover. 

The ordered spectrum is qualitatively similar to the 
one of AF cuprates parent of the high-T c superconduc- 
tors. There the spectrum was interpreted as the superpo- 
sition of two-magnon propagating spin wave excitations 
plus a prominent high energy band due to incoherent 
multimagnon processes, for which no well establish the- 
ory exist. Comparison with CuO suggest that a spinon 
interpretation is reasonable 2 . Furthermore, the relative 
weight of the two propagating spin wave processes re- 
spect to the higher energy excitations decreases as one 
approaches spin liquid states changing the spin or the 
dimensionality. For example in spin-1 2D La2Ni04 the 
spectrum is exhausted by the two-magnon peakH, in 2D 
cuprates 43% of the spectral weight is in the two-magnon 
peakii, here it amounts to only ~ 10% and in S^CuOg 9 
it is zero. This progression gives credit to theories which 
argue that 2D cuprates are close to a spin liquid phase. 
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